Dirac-field model of inflation in Einstein-Cartan theory 
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We present a cosmological model in which a single Dirac field with a flat potential can give rise to 
inflation within the framework of the Einstein-Cartan theory. It is shown that our Dirac-field model 
leads to a nearly scale-invariant spectrum of density fluctuations owing to the spin-interaction which 
naturally arises from the field equations of the Einstein-Cartan theory. 
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I. INTRODUCTION 

Nowadays, it is widely accepted that inflation of the 
early universe [l| is due to one or more slow-rolling scalar 
fields, called inflaton, which typically lead to a nearly 
scale-invariant power spectrum of primordial density fluc- 
tuations with the spectral index n « 1 in good agreement 
with recent, increasingly accurate observations of cosmic 
microwave background anisotropy 2, 3, 4, 5]. However, it 
remains a fundamental task to identify the inflaton with a 
certain particle in promising particle-physics models the- 
oretically and experimentally. Then the question arises: 
cannot other fields drive inflation? Particularly we focus 
on Dirac spinor fields. 

As the scalar field is often necessary in various theories 
of nature, so the Dirac field is essential to a satisfactory 
description of relativistic fermions. Thus, it is of consid- 
erable importance and interest to investigate roles of the 
Dirac field in cosmology. 

On this subject, much remarkable work has been done 
within general relativity by various authors. For exam- 
ple, quantization of a Dirac field and studies of particle 
creation in the expanding Friedmann-Robcrtson- Walker 
(FRW) universe were presented in Refs. [6|, 0,[1|; quan- 
tization of the system where a Dirac field couples to the 
FRW metric was discussed in Refs. d, [l^; the explicit 
solutions of the Dirac equation in the FRW universe were 
derived in Refs. 

Moreover, as for the issue of the accelerating universe, 
it was shown that a self-interacting or nonlinear Dirac 
field can yield negative pressure and thereby accelerate 
the early and the late-time expansion of the universe 
However, when we discuss whether a Dirac field 
can be an alternative to the inflaton, we should pay atten- 
tion to Armendariz-Picon and Greene's work [l5| which 
concluded that, although the de Sitter expansion can be 
driven by a flat potential of a Dirac field, such mod- 
els lead to a scale-dependent spectrum with n — A, and 
thus are incompatible with the observations. As other 
ways not to rely on the inflaton, one can introduce vec- 



tor fields [16|, li2[, or non-Dirac spinors that have some 
unusual properties p^ . 

On the other hand, since Dirac fields have spin- 1/2 
in contrast to scalar fields, it is also of great interest to 
study cosmological effects of the spin. Then, one needs 
gravitational theories that naturally bring spin of matter 
fields to the geometry of spacetime, since in general rela- 
tivity microscopic quantities such as the spin are usually 
neglected. For this purpose, as one of such extended the- 
ories, we will adopt the Einstein-Cartan theory [l9'| in 
which nonvanishing torsion is algebraically equivalent to 
spin of matter fields through the field equation and, as 
a result, a spin-interaction is generated. Generally, one 
can suppose afhne connections in spacetime to be asym- 
metric, thereby naturally obtaining a modified theory of 
general relativity. Since, from the viewpoint of high en- 
ergy physics, general relativity should appropriately be 
extended so as to be able to describe the very early uni- 
verse, the Einstein-Cartan theory can be a reasonable 
framework for considering inflation. We can then expect 
that the spin has a crucial role in inflationary cosmology. 

Many authors have also presented a variety of cosmo- 
logical models based on non-Riemannian gravity [13, Hlj , 
as exhaustively listed in Ref. [25|. The subject of early 
investigations based on the Einstein-Cartan theory was 
mainly to construct exact solutions and to avoid the ini- 
tial singularity in the presence of torsion [l^, [H, [2^ . 
Subsequently, after the advent of inflationary cosmology, 
inflationary solutions were obtained also in the Einstein- 
Cartan models, some of which [2^ [l^, [27l. 28| relied on 



MM- 



'Electronic address: 't omoki@tokai-u.jp| 



spin effects of the so-called spinning fluid 

In this paper, we show that a single Dirac field can give 
rise to inflation within the Einstein-Cartan theory, and 
prove compatibility of the Dirac-field model with the ob- 
servations by calculating the power spectrum of density 
fluctuations of the Dirac field. 

We will take the Dirac Lagrangian as the source of 
metric and torsion because we are only concerned with 
seeking the origin of inflation in the context of particle- 
physics. It should be noted that we will deal with the 
homogeneous Dirac field classically, namely, as a set of 
complex- valued functions that transform according to the 
spinor representation of the Lorentz group and fulfill the 
Dirac equation. We here consider the expectation value 
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of a spinor operator-valued field to be the corresponding 
classical spinor; a possible justification for the existence 
of classical spinors is made in Ref. [l^ . In addition to the 
fact that the classical spinors are mathematically consis- 
tent, as shown in Ref. [15], such classical treatment for 
a physical state can be a good approximation in some 
cases. 

The organization of this paper is as follows: in section 
im we briefly review the description of a Dirac field cou- 
pled to gravity in the Einstein-Cartan theory and derive 
cosmological equations in the flat FRW universe. Within 
this framework, in section [1111 after showing that a flat 
potential of the Dirac fleld leads to the de Sitter expan- 
sion of the universe, we compute the power spectrum of 
density fluctuations of the Dirac fleld. We will observe 
that the obtained spectrum is nearly scale-invariant ow- 
ing to the existence of the spin-interaction. We present 
our conclusion in section IIVI 



II. FORMALISM 

In this section, we summarize how the Einstein-Cartan 
theory describes interactions between a Dirac fleld and 
gravity, and present the cosmological equations that form 
the basis of our Dirac-field model of inflation. 

In order to deal with a Dirac fleld ^ coupled to gravity, 
one requires that the action of the Dirac field be invari- 
ant under local Lorentz transformations, with the help 
of the vierbein formalism (3^ . [33l| in which vierbeins e^" 
and spin connections w^"'' are introduced as the funda- 
mental variables. The vierbein satisfies r]abe^°'ej' = g^i, 
with the Minkowski metric rjah^ where the signature we 
use is ?7ah = diag(-t-l, — 1, — 1, — 1). (The Greek and 
Latin indices denote spacetime and Lorentz indices re- 



spectively.) The spin connection 



is defined by 



where is a covariant derivative 



based on afiine connections of spacetime F''^^ and acts 
on tensors. For the Dirac field, the covariant derivative 
based on the spin connection is defined by 



(1) 



where = («/2)[7°,7^] is the generator of the spinor 
representation of the Lorentz group and the constant 7- 
matrices 7° satisfy the Clifford algebra {7°, 7^} = 277"''; 
the covariant derivative acting on the Dirac adjoint if) = 

-0'I'7O is 



(2) 



In this paper, for explicit calculations, we will choose the 
Dirac representation 



1 

-1 ; ' 



cr™ 

-cr" 



where cr™ are the conventional 2x2 Pauli matrices. For 
later convenience, we also define the additional 7-matrix 
as 75 = i7'^7"'^7^7'^. It is obvious that the covariant 
derivative defined by Eq. ([1]) transforms as a vector un- 
der diffeomorphisms and as a spinor under local Lorentz 
transformations. The coordinate components of 7" are 
defined by 



r 



(4) 



which are shown to give a new set of 7-matrices that 
satisfy the algebra {7'', 7''} = 2^'"''. 

Thus, the Dirac Lagrangian generalized into a curved 
spacetime background is given by 



[i^YD^,^ - {D^lP)-i^'^] - V, 



(5) 



where the term V generically represents a scalar poten- 
tial of the Dirac field including a mass term and self- 
interactions, and consists of arbitrary functions of invari- 
ants generated from -0 and ip; in what follows, we assume 
V = V{s) with s = Tpip for simplicity. 

On the other hand, for gravity, we take the Einstein- 
Hilbert Lagrangian 
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16nG 



(6) 



where R ni, is the curvature of the spacetime given by 



(7) 



In order to take into account effects of the spin of a 
Dirac field on gravity, it is necessary to extend general 
relativity. In this paper, for this purpose, we introduce 
the Einstein-Cartan theory [l9| in which nonvanishing 
torsion C^^^, = 2T^ ^^^^ is related with spin of matter 
fields. (For a very brief review of the Einstein-Cartan 
theory, see Appendix □) 

Within this framework, we consider a system in which 
a Dirac field provides the unique source of metric and tor- 
sion of spacetime. Since the existence of torsion leads to 
the spin-interaction in the curvature ([7|) and the kinetic 
terms of the Lagrangian the conventional Einstein 
equation is modified as 



Gf^p = 8ttG ( T(^,y) 



SttG 



STTGT^r^- (8) 



The additional term —3TTG(t)a4'°'gtii^/'2 with 0" = '0757"V' 
is the spin term prescribed by the Einstein-Cartan theory, 
and can be regarded as a correction due to the spin of 
the Dirac field. Here, the tilde indicates quantities free 
from torsion; G^i, is the Einstein tensor composed of the 
Riemannian (Levi-Civita) connection; ^(^ji/) is the usual, 
symmetric energy-momentum tensor of the Dirac field, 
given by 



T, 



3 



where is the covariant derivative based on the Rie- 

mannian spin connection (1)^°'' = e'''" V^e,y^' , and is 
the torsion-frec Lagrangian of the Dirac field defined by 



- V. 



(10) 



It is interesting that, in Eq. ([8]), the spin correction to 
the energy-momentum tensor of the Dirac field is pro- 
portional to both the metric gfj,i, and the gravitational 
constant G; such a simple form strongly motivates us to 
study the Dirac-gravity system including effects of the 
spin. 

Similarly, the conventional Dirac equation in curved 
spacetime is also modified so as to include the spin- 
interaction term: 



(11) 



where the prime denotes the derivative with respect to 
s. As can be seen from Eqs. ([5]) and ([TT|) . one can arrive 
at the usual general-relativistic equations whenever the 
energy scale of the spin-interaction G(j)a4>°' is negligible in 
comparison with typical energy scales of the kinetic term 
or the potential. 

With the general formalism described above, we are 
now interested to investigate cosmology. Let us consider 
the flat FRW universe, in which the metric is given by 



and the vierbein is chosen to be 



ie;)^diB.g[l,a{t),ait),ait)], 



(12) 



(13) 



where t is the cosmic time and a{t) is the scale factor. 
(See Refs. d, \^ for the nonflat FRW universe.) Then 
we should exclude any spatial dependence of the Dirac 
field for consistency with homogeneity of the spacetime: 
diip = 0. Consequently, we obtain the following Dirac 
equation: 

ij + ^HiP + - 3i7rG0a7"757> = 0, (14) 

where H = a/a is the Hubble parameter and the dot 
denotes the derivative with respect to t. From Eq. 
the anisotropic components Tj^oi) can be shown to vanish, 
which is consistent with isotropy of the spacetime, Goi = 
0. 

Now the energy density and pressure, which must 
be specified to describe cosmological dynamics, can be 
found from the total energy-momentum tensor defined 
in Eq. ^ as 



SttG 



Ptot 



_ (tot) 



sV -V 



SttG 



(15) 
(16) 



It should be noted that (j)a4>°' is always negative, and 
consequently that the spin component in Eq. (|16p con- 
tributes as extra positive pressure. In terms of these 
variables, the cosmological evolution equations are given 
by 



SttG 
AttG 



V - 



37rG 



(17) 



{~2V + isV' - Q'KG(t)a(t)'') , (18) 



which are linked through the conservation law 

Ptot +3iI(ptot + Ptot) = 0. (19) 

One can verify that the conservation law ()19|) is equiv- 
alent to the Dirac equation (fTi|) . Therefore, the system 
that consists of a Dirac field and the flat FRW spacetime 
is self-consistent in the Einstein-Cartan theory as well as 
in general relativity. 

Finally, we define the equation of state w as 

sV -V -iirGdj^ri^ ^20) 



Ptot 
Ptot 



y-37rG(?!)a0"/2 



It is worth mentioning that, by taking the torsion-free 
limit G4)a<i>°' /V 0, we can always reproduce the cos- 
mological equations in general relativity formulated by 
Armendariz-Picon and Greene [lB| • As an important ex- 
ample, the equation of state they found is 



WGR 



(21) 



III. INFLATION 
A. Background 

In this subsection, we show that a Dirac field can cause 
the de Sitter expansion of the background universe on 
the basis of the cosmological equations collected in the 
previous section. 

In general relativity, from Eq. (pij) . one can discuss 
de Sitter infiation driven by a Dirac field, as done in 
Ref. [l3|- In this case, V is assumed to be sufficiently flat 
because wqb. ~ — 1 is guaranteed by 



d\nV 



dins 



« 1, 



(22) 



which is a similar requirement to the slow-roll condition 
for conventional scalar-fleld models of inflation. The con- 
dition ([22| is satisfled if V is asymptotic to a positive con- 
stant for large s. For example, ln(l -I- s"), tanh(ns), and 
s"/(l -I- s)" with positive n are candidates for a potential 
that realizes the de Sitter expansion. 

Also in our case, the flatness of I^ is a basic premise 
for considering inflation. First, in the evolution equation 
for (had)"'. 



dt 



(23) 
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which is derived from Eq. (|14p. the flatness condition 
allows us to ignore the third term on the left hand side 
and hence leads to (j>a't'°' oc a~^. Then it follows that 
all the spin components, which arc incorporated in the 
form of G(j)a4''^ in the background equations presented 
previously, decrease faster than the other terms in those 
equations as the universe expands. 

When the spin components become so small as to be 
negligible, all the background equations in the previous 
section arrive back at those of Ref. [l^ . If the potential 
retains the sufficient flatness until that time, then the 
equation of state eventually becomes w « wgr ~ ^1- 
Thus, a flat potential guarantees de Sitter inflation in our 
model as well. In FIG. [U we have numerically demon- 
strated that w approaches the value of —1 as the spin 
components decrease under a sufficiently flat potential. 

Meanwhile, we can also speculate about the early uni- 
verse before the inflation era. The result that the spin 



components evolve as 



oc a 



also means that these 



become greater in earlier stages of the universe. Thus it 
can be seen that the very early universe is dominated by 
the spin components in our model, and also that, from 
Eq. (|20p . during such an early epoch the Dirac field sat- 
isfies w ~ 1 which corresponds to the equation of state 
for a stiff fluid. Therefore, our model suggests that in- 
flation begins when the spin-dominated era ends, which 
can be seen also from FIG. [TJ Since the Einstein-Cartan 
theory is expected to be relevant to h igh energy physics, 
especially theories of supergravity [s^jISl in which the al- 
gebraic equivalence between torsion and spin is necessary 
for the construction of supersymmetric transformations, 
it is acceptable that the spin components dominate the 
very early universe which should be considered to be at 
high temperature. 

On the other hand, the evolution of s is given by s oc 



during inflation because Eq. p4)) can be recast as 



s + 3Hs + 6iTTG(j)°i:-/5ip = 0, 



(24) 



where the third term on the left hand side arises from the 
spin-interaction, and hence approximately vanishes when 
inflation begins. Consequently, number of e-foldings be- 
tween a time ti and a later time t2 is given by 



iV.lln4^. 

3 S{t2) 



(25) 



Therefore, if one requires w 60 in order to solve the 
horizon problem, then s must change by approximately 
eighty orders of magnitude. In this study we assume the 
potential to be so sufficiently flat as to satisfy Eq. ([^ 
for such a wide range of s regardless of its origin; our 
purpose is to investigate whether the Dirac field with a 
flat potential can be a source of inflation in the Einstein- 
Cartan theory. 

Similarly to the scalar s, the pseudoscalar 55 defined by 
S5 = V'75'0 evolves according to S5 cx during inflation 
because S5 fulfills the evolution equation 



where both the last and the penultimate terms on the 
left hand side are negligible under our assumption. We 
also note that, in terms of the rescaled field defined by 
"if = a'^^^ip, these scalars, s and 55, can be rewritten as 



'\l/75^' respectively, where and 



^'75^' are constant as long as the sufficient fiatness of the 
potential holds. 




S5 + 3Hs5 + 2iV'(f>" + QinGctPs = 0, 



(26) 



time 



FIG. 1: Typical time evolution of the equation of state w. 
The solid curve represents w for the solution of the cosmo- 
logical equations (I14|l and (llSp with the potential of the form 
ln(l + 8^), and shows that its value which is initially w ~ 1 
comes close to —1 as time progresses. The dashed curve repre- 
sents the simultaneous evolution of the spin components in w, 
(-37rG<^a(j!>''/2)(l/ - 37rG(^a0''/2)"\ disappearing with time. 



B. Perturbation 

In this subsection, we discuss perturbations of the 
Dirac field that brings about de Sitter inflation via the 
dynamics explained above. 

The conventional scalar-field models of inflation [l| 
have an important feature of predicting a nearly scale- 
invariant spectrum of density fluctuations in excellent 
agreement with recent observations 0, S 0, 0] • Our next 
task is thus to examine the consistency of the Dirac-field 
model of inflation with the observations, namely, whether 
our model can derive a nearly scale-invariant spectrum, 
by computing the power spectrum of density perturba- 
tions of the Dirac field. Although a proper analysis for 
this purpose should be based on gauge-invariant pertur- 
bation theories [H, [13, [H, where both spacetime and 
matter fields are perturbed, we will not consider the met- 
ric perturbations but only the perturbed field 5^ for sim- 
plicity. 

Our key strategy is to perturb the spin components of 
the background equations. Whereas Armendariz-Picon 
and Greene [l^ showed that, within the framework of 
general relativity, Dirac-field models of inflation with a 
flat potential lead to a scale-dependent spectrum and 
hence are inconsistent with the observations, our model 
based on the Einstein-Cartan theory possesses in the first 
place the spin-interaction, which opens the possibility of 
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improving their result. It should be noted that the terms 
related to the spin in our model appear as a consequence 
of a natural extension of general relativity. As discussed 
previously, the inflationary expansion itself can be shown 
to occur by simply assuming the flat potential, whether 
the Dirac field has the spin-interaction or not; however, if 
the spin terms exist, then their fluctuations must also ex- 
ist. We will show below that our model predicts a nearly 
scale-invariant spectrum owing to these fluctuations. 

In what follows, we mark the background quantities 
with B, and for convenience, use the identity 

cp^r = -{i'^f - {th5^f = -s' + si (27) 

The scalar s = ipip and the pseudoscalar 55 = ip^^ip are 
now perturbed as 

s = ^pBipB + i'B^i' + Si^tpB ^ SB + 5s, (28) 

55 = i>Bl5i>B + i>Bl5SlJj + SlJjj5i>B = sf -f- (5S5.(29) 

The equation of motion for Sip is then given by 

^3^G [{sB - sf 75)(5V' + (Ss - <5s575)V's] = 0, (30) 

where m = V' and we have ignored V" because of the 
flatness of V, from which it follows that m is approxi- 
mately constant. Again we note that the term propor- 
tional to G is the correction due to torsion or equiv- 
alently spin; i.e., one can take the limit G — > as a 
method to render the perturbed equations torsion-free in 
the dynamics of perturbations, while, in the background 
equations, the torsion-free equations are reproduced as 
a result of rapid disappearance of the spin components 
of the form G(j)a4'"'- On the other hand, taking account 
of the fact that general relativity provides a successful 
description of the present universe, one can consider the 
correction coming from torsion to be generally small. In 
other words, we can regard the term proportional to G 
as a first-order correction to general relativity, and will 
ignore higher-order corrections in G whenever the quan- 
tities derived below include such terms. 

In the limit G ^ 0, the torsion-free equation of motion 
is found to be 



ij'^Df^Stp — mStp — 0, 



(31) 



which can be solved analytically For H w const., 

the plane- wave solution Us(k,t) of Eq. (|?T|) is generally 
given by 



r(i.2)/ 



(32) 

where rj is the conformal time and H^''^' {—kr]) are the 
Hankel functions of the first and the second kind with 
u = 1/2 — im/H « const.; and (3^ are arbitrary 



constant two-spinors. The solution of Eq. pTjl can then 
be expanded as 



Si; = 



d^k 
(2^ 



[us(k,t)as{'k)e 



+f,(k,t)6t(k)e- 



(33) 



where the mode function Vs (k, t) corresponding to nega- 
tive energy solutions is of the same form as Eq. ([32)1 . 

Now, for the purpose of investigating the changes due 
to the spin from the torsion-free case, it is useful to intro- 
duce scalar functions A(k, t) and B(k, t), and to express 
the plane-wave expansion of the solution of Eq. (j30p in 
terms of the analytical solution (|32p as 



Sip = 



Sk 



^ ' S—1 

fi3(k,i)i,,(k,i)6t(k> 



— ik-: 



(34) 



in which expression, A and B play a role to represent 
the modifications to the torsion- free case; i.e., by sub- 
stituting Eq. p4|) into Eq. pO|) . one can see that the 
evolution of A and B is governed by the spin-interaction 
part of Eq. ([50]) . and also that both A and B are con- 
stant in the limit G — > 0. Here, and 6s are the particle 
and antiparticle operators satisfying {as(k), a^, (k')} = 
{6s(k), 6^, (k')} = (5ss"5'^(k — k') under an appropriate or- 
thonormality. Then, the vacuum expectation value of the 
square of (5-0 is 



{STpSip) = 



(Pk 
(27r)3 



|Sp5^S,(k)«,(k). (35) 



Next, we characterize density perturbations 5p by the 
variable 



5p ^ 



(36) 



This quantity is the one employed in Ref. [l5|, de- 
fined analogously to the gauge-invariant Bardeen vari- 
able, which is conserved for adiabatic perturbations on 
sufficiently large scales [S^H^]- For our model, we have 
p + p = msB and 



Sp = mSs + 3t:G{sbSs — s^ Ss^), 



(37) 



from which we compute the power spectrum V of the 
variable ^ on comoving scales 1/fc, defined by 



(C(x,t)C(x + r,i)) 



dk sin kr 
k kr 



V{k). (38) 



We have supposed that the Fourier transform of the cor- 
relation function on the left hand side depends only on k 
and not on k itself, for consistency with the isotropy of 
the background spacetime. 

In order to evaluate the isotropic power spectrum of 
(C^), following the method of Ref. IS], we apply the Fierz 
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transformation to (C^) and expand the vacuum expecta- 
tion values of the quadratic term of 6s and Ss^ in terms 
of perturbation bihnears. For example, the vacuum ex- 
pectation value of Ss^ can be expressed as 



(39) 



After that, as noted in Ref. [l5|, we concentrate on 
the perturbations of the scalar bilinear form (di/jSilj) 
and discard the remaining vectorial bilinears that may 
break isotropy of the power spectrum because we are 
only concerned with determining the amplitude and k- 
dependence of P. In this way, V is obtained as 



r{k) 



\B\^k^ 



J2^sVs, (40) 



where Vg is the rescaled mode function defined by Vg = 
a^/'^Vs and Ci ^ 2^ b - (^b75^'b)^/*b*s is a con- 
stant. We have here ignored the terms of 0{G^) as men- 
tioned before. In comparison with the result of Ref. p^ . 
the extra factor |_Bp appears in addition to the first-order 
correction in the parentheses in the above equation. The 
effective part of B that contributes to the isotropic power 
spectrum fulfills the following equation: 



B 
B 



'SinG 



J2 (si3^s7°K 



+i>BV,Zl°i'B - ^BToVsZl°l5i'B), (41) 



which is obtained by multiplying Eq. (|30p by Sip with 
Eq. (j34p and evaluating in a vacuum. 

Now we have to determine the mode function Vs and 
B. We here take the Bunch-Davies vacuum as an appro- 
priate initial vacuum state [33| , and choose the constants 
in Vs so that Vs coincides with the Minkowski solution 
Vs cx e''^'' on sufficiently short scales k 



oo: 



Vs = 



2 I e"""'^^" Hll\~kf])Lu' 



(42) 



where is the constant two-spinor defined by to^ = ^a- 
With this choice, in order to solve Eq. (|4T|l . we use the 
explicit background solutions during inflation, 



a = a*e 



(43) 
(44) 



where (po and xo £^re arbitrary constants. Moreover, for- 
tunately, we can make use of the fact that VsVs is ap- 
proximately constant with respect to z = —krj. Roughly 
speaking, the relation \d{Y. ^ ZVs) / dz\ < lE.^s^sl 
holds for z > 1 and near 2 = independently of to. 
In fact, the asymptotic formulae of the Hankel functions 

HH jead to J2s ^sVs 2{m/H)z-^ for large z and 

^s^sVs ^ —2tanh{7Tm/ H) for small z. We can safely 
employ this relation because, as discussed below, we only 
need to find the asymptotic behavior of on large and 
small scales. Then Eq. (PT|) can be integrated to give 



\B\ 



2Y.^VsVsk^ i \a,H 
7rTO„ f 3 



, -.21771/ H \ 3 



cosh 



H 



- i sinh ^ + z2''2F3 + l,iy+^;-z^ 

^-4i7n/H / l\ / \\ / I 1 3 

-— . + - r 2.+ -U^3U+2,2.-H-;.,2.,2.+ -;-z^ 



1 



- h.c. + C2 



(45) 



where 2^3(01, 02; foi, ^2, ^3; x) is the regularized hyperge- expansion of 2-^3(01, 02; 61, 62, ^3; x) for |a;| 00, 
ometric function defined by 



2i^3(ai,a2;&i,&2,63;a;) 



2i^3(ai,a2;&i,^'2,&3;a;) 



with the Pochhammer symbol {a)k = T(a + k)/T{a), and -| 

we have dropped the higher-order terms in G. The inte- ^(^i ^ 01)^(^2 ^ ai)r(fe3 - ai)r(a2) 

gration constant G2 should be chosen so that 1 r(ai — 02) (—a:)^"^ 

on sufficiently short scales k ^ 00. Using the asymptotic ~'~r(6i — a2)r(62 — a2)r(fe3 — a2)r(ai) 



r(a2-ai)(-x)-'^i 

si)r(a2) 

(47) 
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with c = (fli + a2 — bi — b2 — b^i + l/2)/2, we have 



C2 



sech- 



27rm 



<y5oXor(-i') 



r(--i)r(i 



2im/H 



2v) \a^H 



.(48) 



The power spectrum T' outside horizon is obtained 



by evaluating Eq. pOl) on large scales ^ with 
Eq. (|45p . in which we can use the asymptotic formula 

2^^3(01, a2; b^MM; x) ^ (r(6i)r(62)r(63))"' for x -> o. 

Since the observed spectral index at horizon crossing is 
defined by n — 1 = d\nV / d\nk\k=aH , wc subsequently 
estimate 7^ at fc = aH. For simplicity, let us here uti- 
lize an approximation m/ H <^ 1, which is allowed by the 
flatness of V . Thus we find 



J 



V\k=a 



H 



^(2C. + -(^U + xS^o) (^(37 + ln8-l) 



J- A(3^ + ln8-i; 
4 37r 



In 



-|-(37 + ln8 - 2) + ^(27 + ln4 - 3)] 



(49) 



r 



where 7 is Euler's constant, 7 
spectral index is 



0.577. Therefore, the 



m 



iivlxo-xl-fo) (7^(37 + In 8- 1) - 



97r 



1 
2 
(50) 

which is our final result. It should be emphasized that 
the fc-dependence of V is naturally suppressed by the 
condition in/ H <^ 1 which follows from the flatness of 
the potential, and hence that n is nearly equal to 1. We 
also note that the running index a = dn/ dink is further 
suppressed, of the order O(to^). Therefore, we are led 
to the conclusion that our Dirac-field model of inflation 
can predict a nearly scale-invariant spectrum of density 
fluctuations in agreement with the observations. 



IV. SUMMARY 

We have shown that the Dirac-field model of inflation 
leads to a nearly scale- invariant spectral index consistent 
with the observations, by naturally extending the theo- 
retical framework beyond general relativity. It is usually 
believed that general relativity does not hold in the very 
early universe and then should be extended appropri- 
ately to high energy physics; the Einstein-Cartan theory 
adopted here is one of such extended theories. 

In the framework of general relativity, Armendariz- 
Picon and Greene !15| found that a Dirac field with a flat 
potential can give rise to the de Sitter expansion of the 
universe, but concluded that the Dirac field itself cannot 
be an alternative to the conventional inflaton field as the 
unique source of inflation because the spectral index ob- 
tained from density fluctuations of the Dirac fleld, n — 4, 
is in strong disagreement with the observations. The key 



ingredient for improving their result, namely, obtaining 
n « 1, is to introduce a spin- interaction, which naturally 
appears in the Einstein-Cartan theory, into the dynamics 
of the Dirac fleld that has the inflationary potential. 

Because of the existence of the spin-interaction, the 
new terms of the form G(j)a4>°' rnust be added to both 
the cosmological background equations and the equation 
of motion for the perturbed fleld presented in Ref. ^15| . 
However, we have seen that, according to Eq. ((23)) . the 
spin terms in the background equations decrease fast as 
the universe expands; therefore, during inflation, we can 
ignore the effects due to the spin on the background. 
Without any change in the basic idea that a flat potential 
of a Dirac field leads to an inflationary expansion, we 
have been able to gain the possibility to solve the problem 
of the spectral index, which is an issue in perturbation 
theories, not in background dynamics. 

The spectral index obtained in our model, n = 1 + 
O(m^), is nearly scale-invariant by virtue of the flatness 
of the potential, which is similar to the situation in the 
conventional infiaton models. In this regard, so far, the 
Dirac-field model does not provide such novel features as 
the typical inflaton models do not have. However, it is 
important to recognize that a Dirac field can drive infla- 
tion of the universe. As is well known, the spinor fields 
are indispensable not only in the description of relativis- 
tic quantum fields, but also in the context of supersym- 
metric unification of all fundamental interactions at high 
energy scales. Therefore, in the construction of realistic 
cosmological models containing various matter fields and 
interactions between them, attention should be paid to 
the behavior or properties of the spinor fields that signifi- 
cantly affect the geometry of spacetime and consequently 
can have a central role in the evolution of the universe. 
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APPENDIX: EINSTEIN-CARTAN THEORY 

For completeness, in this appendix we briefly review 
the Einstein-Cartan theory [1^. The Einstein-Cartan 
theory is a natural extension of Einstein's gravity theory, 
and is one of theories that give a dynamical role to both 
spin and mass of matter. 

The spacetime in this theory is described by Riemann- 
Cartan geometry known as a generalization of Riemann 
geometry to include torsion. In the Riemann-Cartan ge- 
ometry, from asymmetric afhne connections T^^^^, the 
covariant derivative for tensors is defined by \7uV^ — 
d^V^ + V^p^V. The connection is constrained by the 
metricity condition V pg^i, = 0, which is postulated in 
order for a local Minkowski structure to be guaranteed. 
The curvature tensor is constructed from such connec- 
tions as 

R''„,.^d,Tf^^-d.T'^,^+T\^T\^~T\^T\^. (A.l) 

The difference from the Riemann geometry is that V 
is asymmetric. Actually, if we demand the connection be 
symmetric, it can be fixed as the well-known Riemannian 
connection 

= Ig'^d^g,^ + d^g„, - d^g^,). (A.2) 

There is, however, no a priori reason that we suppose 
V to be symmetric in general. In the Riemann-Cartan 
geometry, the antisymmetric part is kept as 

Cpiu = 2r''[^yj = r''^j, — Vjj^, (A. 3) 

which can be shown to transform as a tensor, a purely ge- 
ometrical quantity. Since the infinitesimal parallelograms 
do not close in this spacetime and the closure failure is 
proportional to C'^^j,, this tensor serves as the torsion of 
the spacetime. 

These geometrical quantities, the curvature and the 
torsion, can be understood from the local Poincare gauge 
theory [H, [U HI, HI, HI] in which vicrbcins e''^ and spin 
connections uj^""^ are introduced as the gauge fields of the 
theory. In terms of the gauge fields, we can define the 
translational field strength corresponding to the torsion 
as well as the rotational field strength corresponding to 
the curvature [isj . 

From the correspondence between a coordinate basis 
and a tetrad, the vierbein satisfies 

S.e^'^ + T^p.eP'' + c^/^e^, = 0. (A.4) 

The above equation guarantees a conversion V ^V^ — 
e'^g^D fj^V"" , where Z?^ is the Lorentz covariant derivative 
based on w, 

By virtue of the metricity condition, T^^^ can be 
decomposed into the Riemannian piece and the non- 
Riemannian piece as 



from which it is clear that the connection reduces to the 
Riemannian one when torsion vanishes. Such decompo- 
sition holds also for the spin connection with the aid of 
Eq. jXl. 

The basic field equations of the Einstein-Cartan theory 
are derived from the simplest Lagrangian of a gravity- 
matter system 

£(e"„<^<^,9,,^) = --^+C,r.ie^,,^,D^ip\ (A.6) 

where R = e'^^e'^jji?"''^^ is the scalar curvature in the 
Riemann-Cartan spacetime and ip generically represents 
matter fields minimally coupled to gravity. Here, by mak- 
ing use of the relation l|A.4|) . not only the matter La- 
grangian Cm, but also the Einstcin-Hilbert Lagrangian 
can be expressed as a function of e^,^ and uj^"'''- 

For a derivation of the field equations, it is useful to 
adopt the so-called Palatini approach where e^^ and uj^"''' 
are treated as independent variables. Varying the La- 
grangian (|A.6p with respect to both, one obtains 

^ V - ^ V, (A.7) 



= -I&ttGS^,,, (A.8) 

where = e'^j^R^'^^^ is the asymmetric Ricci tensor and 
Sab is the generator of the Lorentz group. (The capital 
indices are spacetime or spinor indices.) 

One of the resulting equations, Eq. (jA.7p . is a gen- 
eralized version of the familiar Einstein equation in the 
sense that the Einstein tensor G,," which consists of the 
metric-compatible and asymmetric connection is related 
with the canonical energy- momentum tensor T^°. The 
other equation (jA.8[) exhibits that the spin density 5''^^ 
of matter induce torsion of the spacetime. 

Since the second field equation (jA.8|) is algebraic, we 
can substitute everywhere spin for torsion. Then, accord- 
ing to Eq. (jA.Sp , the non-Riemannian part of the covari- 
ant derivatives in Cm produces an interaction between 
^^ab 'fi^ called a spin-interaction, in the equation of 
motion for ip. 

Moreover, Eq. (jA.5|) can also be applied to splitting 
the Einstein tensor G^'^ defined by Eq. (|A.7p into the 
Riemannian piece and the non-Riemannian piece. There- 
fore, all the terms including the torsion in Eq. (jA.Tp can 
be interpreted as a spin correction to the usual, symmet- 
ric energy-momentum tensor ^(^1^) through the following 
expression: 

G^, =87rG(f(^,)-l-T^r")), (A.9) 

where G^^ is the Einstein tensor composed of the Rie- 
mannian connection f ''^j^. The explicit form of T^it'""^ is 
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given by 



(A.10) 



where T^l^™'' arises from the couphng of </3 to S^^^^ in the 
covaxiantized kinetic terms of Cm.. 



The Einstcin-Cartan theory is constructed in this way. 
Since, as can be seen from the above equation, the spin 
squares contribute to the energy-momentum in the form 
proportional to G, the predictions of the Einstein-Cartan 
theory deviate from those of general relativity only when 
the matter fields are coupled to gravity at high energy 
scales. Therefore, it is reasonable to take the Einstein- 
Cartan theory as a suitable framework for considering 
early stages of the universe, in which general relativity is 
usually not believed to be valid. 
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